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Fluctuation Theorems are statements about the entropy of systems far from thermal equilibrium. 
In this Letter relativistic Fluctuation Theorems for Brownian motion are presented and proven. 
Though there is a known discretization dilemma leading to a certain class of Brownian processes, 
the Fluctuation Theorems are shown to cover the entire class, and to determine the physical process 
within this class. In the general-relativistic case the heat responsible for entropy production is the 
exchanged 4-momentum projected on the local time axis of the heat bath. 
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Introduction. — The physical basis of the direction of 
time is being discussed since Boltzmann's H-Theorem in 
1872. A priori, this thermodynamic arrow of time has 
to be distinguished from the possibility of a prime di- 
rection of time defined by the monotonous expansion of 
our Universe 0. Though we know by now that due to 
its dominating Dark Energy component our Universe is 
very likely to expand forever, a direct connection between 
these arrows of time seems unlikely 0, 0, , leaving the 
fascinating cosmological arrow of time in a speculative 
stage 5, 6]. Highly advanced is the situation of the arrow 
of time in thermodynamics. Since 1993 Fluctuation The- 
orems (FTsjhave been derived for very general classes of 
systems II II H El El E3 Not only are these FTs 
capable of proving the second law of thermodynamics for 
the entropy AS produced by an ensemble, 

(1) 



AS > at any time, 



they even allow statements about the probability of "vi- 
olations" of QJ by a single member of the ensemble in 
finite regions of space and time. For the steady state of 
sufficiently chaotic dynamical systems, FTs of the form 



P(As = +a) _ a 
P(As = -a) " eXP fcr? 



(2) 



valid for any a have been proven |7J,|a,|9j. Here P(As = a) 
is the probability to observe a production of entropy 
equal to a along a trajectory segment. In the context 
of dynamical systems, the dimensionless entropy As/fee 
is defined by the contraction of phase space. The Steady 
State FT (0) was derived in [l(j for a stochastic sys- 
tems in contact with a bath at temperature T. For this 
stochastic systems the entropy As = AQ/T is caused by 
an energy transfer AQ into the bath. In 0, 0, 0] a 
certain entropy production As is averaged exponentially 
to find Integral FTs of the form 



-As/k-t 



1, 



(3) 



which remain valid even for non-stationary states in the 
presence of time-dependent external forces. Angle brack- 



ets denote averaging over trajectories by path integration 
H3- Equation © implies ©, and Eq. © implies (QJ by 
virtue of the Jensen inequality using AS — (As) . 

Recently, the 1905 publications of Einstein on Brown- 
ian Motion Q and Special Relativity E3 were unified 
flil UtI ITsf . Based on that, the results of this Letter are 
threefold. First, for the relativistic Brownian processes 
E3 and H3, we reconcile the FTs 101 and Jjjj, which 
have become a paradigm of nonequilibrium physics, with 
Special Relativity. FTs for a similar process put forward 
in El] can be derived readily by an analogous thread of 
reasoning. 

In [13 and E3 it was pointed out that the relativis- 
tic time dilation leads to multiplicative coupling, neces- 
sitating a careful choice of the three discretization rules 
(according to Ito, Fisk-Stratonovich, as well as Hanggi- 
Klimontovich) , since they lead to physically different pro- 
cesses. As second central result we show, that there is 
one relativistic Steady State FT (J2J) and one relativistic 
Integral FT J3J valid for all choices. Furthermore, we 
shall find the physically correct expression for the en- 
tropy production following from relativistic FTs when 
the Hanggi-Klimontovich discretization rule is applied. 

The third result of this Letter is a general-relativistic 
Integral FT for the Cosmological Standard Model. This 
exposes clearly the role of cosmic expansion in entropy 
production. We shall identify the entropy production 
which is solely due to the Hubble expansion of space. 
Such entropy producing processes dominate when the ex- 
pansion rate of the Universe exceeds the particle scatter- 
ing rate, for instance in an early inflationary phase after 
the Big Bang. 

Relativistic Brownian Motion. — To avoid technical- 
ities we consider first the one-dimensional motion of a 
test particle in a heat bath. The generalization to higher 
spatial dimensions is straightforward. In the high tem- 
perature limit, the mean squared velocity of this particle 
with rest mass m may no longer obey the nonrelativistic 
law v 2 = k^Tjm since the finite speed of light defines an 
insurmountable upper bound. The authors of [lfl . lrRll^ | 



2 



have set forth the relativistic Brownian motion giving rise 
to relativistic velocity distributions in both, the language 
of stochastic differential equations (relativistic Langevin 
equations) and the language of probability densities (rel- 
ativistic Fokker-Planck equations). Following \vH Hi| . 
the generalized deterministic force Fa in the rest frame 
of the heat bath is 



dpd — F<j dt = —vp dt , 



(4) 



so that the time scale of dissipation is 1/v. Equation (0J 
is of the familiar form |20j | with the nonrelativisic momen- 
tum mv substituted by p = p 1 = mv/ 'yl — v 2 /c 2 , which 
is the spatial component of the relativistic momentum 
vector p a . As is common, Greek indices refer to tempo- 
ral (a = 0) and spatial components. The signature of 
the Minkowski metric tensor is r] a p = rf 1 ^ = diag(— 1, 1). 
Moreover, Einstein's summation convention is invoked 
throughout. Since the rest mass is not changed in col- 
lisions, p a p a = — (mc) 2 = const, the change in the mo- 
mentum vector dp a is always "orthogonal" to p a in the 
sense of 



p a dp° 







(5) 



This means that the classical particle cannot leave its 
mass shell p a p a = —(mc) 2 , which is nothing but its dis- 
persion relation, 



E = p°c=^(mc 2 ) 2 + (pc) 2 . (6) 
The general solution of (0 is the projection 

of an arbitrary Lorentz vector It is readily seen that 
Eq. Q , valid in the rest frame of the bath, implies 



dp2 = —mis 



P Pf3 

(mc) 2 



7 bath 



dr 



(7) 



with the bath velocity vector f(, ath and the particle's 
proper time r. Equation Q is the generalized Lorentz- 
invariant deterministic part of the Brownian motion |23| . 

The description of relativistic Brownian motion is com- 
pleted by Lorentz- invariant stochastic changes dp" of the 
momentum caused by the impacts of the surroun ding 
heat bath at temperature T. The derivation in 0, Hsf 
rests on two ideas: The relativistic momentum is the 
proper quantity performing a Wiener process, since it is 
physically exchanged and additive, whereas the velocity 
is well-known to be not additive in Special Relativity. 
The second postulate demands that the distribution is 
a Gaussian in the instantaneous rest frame of the par- 
ticle. This connects the relativistic Brownian motion to 
the nonrelativistic case. These principles determine the 



exchanged momenta dp" do be distributed according to 
(cf. Eq. (35c) in [13) 



P co n(p^dPs 



0~ (Pfidtf) dp«dp s 
■ exp - 



40dr 



(8) 



2Vn¥dr 

The Dirac distribution 5(ppdpP) guarantees that the 
mass-shell condition (01 is also fulfilled by the stochas- 
tic impacts, since they are elastic. While the relativistic 
momentum p is additive and unbounded, the velocity is 
restricted to the open interval (— c, +c). This can be seen 
by the elegant relation v/c 2 = p/E in the rest frame of 
the bath, which is equivalent to 

pc 



dx = 



dt . 



- p- 



(9) 



The bath temperature T is defined by the Einstein rela- 
tion for the momentum diffusion constant S> (cf. Eq. (59) 
in |l7|1: 



k^Tmv 



(10) 



The Einstein relation is known to apply far from equi- 
librium 0| . This is not surprising since (|10|l relates the 
bath temperature to the strength of stochastic impacts 
S>, and the particle damping rate v. These quantities de- 
fine the coupling between the bath and a single particle, 
be it considered as part of an equilibrium ensemble or 
not. 

Relativistic Fluctuation Theorem. — We have now the 
manifestly Lorentz-invariant Langevin equation 



dp a = dp% + dp" 



(11) 



with the deterministic part given by J7J and the stochas- 
tic part described by (jHJ) at hand. Specifying to the 
rest frame of the bath, vg ath = (c, 0), yields 



dp = —vp dt + dp s 



(12) 



The probability density of the exchanged momenta dp s 
is found by integrating out the dpg-component in JHJ, 
cf. 0: 



-P co n(p,dp s ) 



exp 



-dp s 2 /4^ v /l 



(mc) 



?dt 



2n/^V / 1+(^ 



(13) 



The Eqs. iJHJl, lfT2"|l and ifLljl establish the relativistic 
stochastic motion of the Brownian particle in phase 
space. The corresponding transition probability is eas- 
ily determined by the discretization rule, for example in 
the case of Ito: 



Pu 



x i ► x ~\~ dx 
p i— > p + dp 



5 [dx 



pc dt 



(mc) 2 +p' 



2V^d^/l + I ^ 



(dp + vpdt) 
xexp ; - ■ (14) 



40, /l 



(mc) : 



r&t 
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It has been clarified in |13( that entropy fluc- 
tuations are caused by the particle entropy 
s p = — fc B lnP(x(t),p(t),i) with the particle's nonequi- 
librium phase space density P(x,p,i), and the entropy 
s m of the surrounding medium at temperature T. From 
the equation of motion (cf. Eq. (7) in [l3|) for s p , 

ds p = ds p |A=o + Afc B dln£J , (15) 

we have isolated the second term which depends on the 
discretization rule applied to 112[l : Hanggi-Klimontovich, 
Fisk-Stratonovich, or Ito correspond to A = 0, |, or 1 
respectively. 

The entropy production ds m in the bath follows by 
contrasting the probability of a trajectory (x,p) t f with 
its time-reverse = (x(tf — t),—p(tf — £)) to 

extract the time-asymmetric part causing dissipation 
(cf. 0| and references therein): 

ds m = fcB In t , rans = —— — Afce dlnP . (16) 

PLns 

From the Eqs. I|15|) and URifl we see that though the 
relativistic motion of the Brownian particle is physically 
inequivalent depending on A, the fluctuations of the total 
entropy s — s p + s m are independent of A. The Eqs. i|15|) 
and l |l 6|l can readily be integrated by the method de- 
veloped in over a finite time interval to the Steady 
State FT J2J an d for time-depend states to the Integral 
FT Therewith we have proven relativistic FTs that 
are unaffected by the discretization dilemma. 

We are now in the position to address the physical 
choice of A by virtue of the FT. In the nonrelativistic 
regime, E is dominated by the rest mass so that the 
second term in (|16|) vanishes for all A. At arbitrary rela- 
tivistic energies © the Hanggi-Klimontovich rule, A = 0, 
yields the correct expression for the entropy 

ds m = ~ , (17) 

which is produced in the heat bath. 

Generalizations in the framework of Special Relativity. 
- To generalize the FTs to n spatial dimensions, mo- 
mentum and force in Eq. J3J are simply substituted by 
their spatial vectors and the Greek indices in the Lorentz- 
invariant Eqs. J7J and JSJ take values up to n. After inte- 
grating out the temporal component p , the distribution 
IjlHJI is found to contain a quadratic form A instead of 
the square in the exponent (cf. Eq. (15) in 1 S| ) with 
tensor components 

c 2 

A H = S i] - -jpPiPj ■ ( 18 ) 

The FTs follow using the obvious fact that p is eigenvec- 
tor of A. 

No complications are caused by allowing an inhomoge- 
neous heat bath, where the dissipation rate v is a function 




FIG. 1: A sketch of spacetime showing a spatial slice of the 
heat bath at fixed time and the world line of a Brownian 
particle in a (locally) expanding Universe. 



of space and time. The dissipation rate v may also be an 
even function of the momentum, ^(|f>|). This of impor- 
tance since we may not expect the embedding medium 
to behave as a Newtonian fluid at relativistic energies. A 
temperature varying in space and time does not pose a 
problem, nor does an additional time-dependent external 
force. In this case the produced heat dQ = Tds m is the 
loss of particle energy in collisions, — dE reduced by the 
work extracted via the external force. 

Generalizations in the framework of General Relativ- 
ity. — The monotonic increase of entropy is a funda- 
mental principle of physics and the Universe is known to 
expand, as was discovered by E. Hubble in 1929. The 
discussion whether there is a direct connection between 
these observations has newer stopped 0, IH 0, H, E| 
Therefore we aspire a formulation of the FT consistent 
with General Relativity, but we restrict ourselves to the 
class of Friedmann-Lemaitre models, which describe a 
spatially homogenous and isotropic, expanding or con- 
tracting Universe. The corresponding line element (as 
given by the Robertson- Walker metric) is —di 2 + dr 2 . 
The important difference compared to Special Relativity 
is that the spatial part, dr 2 , is scaled by a time depen- 
dent factor R(t) describing the expansion or contraction 
of the Universe: 

dr 2 =R 2 (t)hij{Ode d& . (19) 

The Latin indices describe spatial components numbered 
by 1 to n. We do not have to deal with the details of 
the metric tensor h describing the spatial geometry. The 
result will be valid for arbitrary geometries. The expan- 
sion rate H(t) = R(t)/R(t), named Hubble function, is 
one of the most important quantities in cosmology and 
its present value is a direct observable |2l|. The typical 
frame for a cosmic heat bath is the frame of the cosmic 
microwave background. 

In General Relativity, the correct equations of motion 
include the covariant differential Hp of the momentum. 
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(Denoting by p the 4-vector, the components of Dp are 
Dp a = dp a + T" p^dx u .) Its spatial components replace 
the left hand side of (|12|) and can be split up into a spa- 
tially covariant part, WDp, and a contribution due to 
the time-dependent scaling: 



Dp = (n) Dp + #(t)p . 



(20) 



Therefore the covariant Langevin equation, generalizing 
Eq. lO to be valid in an expanding or contracting Uni- 
verse of arbitrary spatial geometry, reads 



(n) 



Dp = MIIpIU) + #(*)] p dt + ( " )d p 



(21) 



The distribution of the stochastic impacts (™'Dp s is found 
after substituting the Euclidean metric Sij by hij in (I18|l . 
Applying the time-reversal map we find that Eq. ltT7|) 
gains a second term: 



ds r 



dE 
T 



|PC|| 



dlni? 



ds 



static 



ET 
(P,dr) 



H 



T 



(22) 



The numerator (p, dr) in H22|) is the canonical line inte- 
gral (canonical one-form) in phase space. The Integral 
FT J2J extends to an expanding (H > 0) or contracting 
(H < 0) spacetime when this second term is taken into 
account. It has a clear geometric interpretation: The 
Hubble function is the external curvature of space, 



DN = H dr 



(23) 



with ./V being the time-like normal vector to the space 
of the heat bath as depicted in Fig. This permits the 
second term in Q22JI to be written as 

T 

Since the particle energy E = p° = —po = — (p, N) is the 
zero component of the 4-vector p, the first term in 122|) 
equals the differential 

d(p, N) _ (Dp, N) + (p, DN) 
T T 

such that the sum of both terms is the projection 



ds„ 



(Dp, AO 
T 



(24) 



As described by Eq. \2'dty . the geodesic flow N may be 
convergent or divergent, corresponding to a contracting 
or expanding Universe, but this does not imply a change 
of sign for the entropy production ds m . 

Conclusions. — Relativistic FTs have been established 
that remain valid in the regime of high temperatures or 
low masses, mc 2 <C k^T. By investigating the entropy 
production for particle and environment separately, we 



could determine the physically correct discretization rule, 
which had not been possible sofar by relativistic invari- 
ance alone. These FTs were found to extend in the frame- 
work of General Relativity. Such a formulation reveals 
the influence of a time-depended gravitation field on the 
local entropy production. The dissipated heat Tds m is 
the exchanged 4-momentum Dp projected on the local 
time direction of the heat bath. 

On the theoretical road ahead, one may expect Inte- 
gral FTs to hold for arbitrary time-dependent and inho- 
mogeneous fields, such as gravitational waves, when the 
concise expression H24[l is applied. For the process orig- 
inally introduced in 1161. the weaker inequality has 
been proven recently |22( under general conditions. Ex- 
perimentally, the relativistic FT is not only subject of 
high energy physics and cosmology. An ultra relativistic 
FT can be tested with a high-precision spectroscopy ex- 
periment by shining with a laser on an exited granulate 
of glass beads, so that the granulate serves as a heat bath 
and the photons are the relativistic "Brownian" particles. 

I am grateful for comments of J. Dunkel, S. Herming- 
haus, U. Seifert and M. Brinkmann . 
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